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ABSTRACT 

A reliable method is presented for calculating the flontfield about a 
cascade of aroitrary 2-0 airfoils. The method approximates the 
three-dimensional flow in a turbomachinery blade row by correcting for 
streamtube convergence and raoius change in the throughflow direction. The 
method is a fully conservative solution of the full potential equation 
incorporating the finite volume technique on a body-fitted periodic mesh, 
with an artificial density imposed in the transonic region to insure 
stability and the capture of shock waves. Comparison of results for several 
supercritical blades shows good agreement with their hodograph solutions. 
Other calculations for these profiles as well as standard NACA blade 
sections indicate that this is a useful scheme for analyzing both the design 
and off-oesign performance of turbomachinery blading. 

INTRODUCTION 

Some very useful tools for the study of transonic flows have been 
developed during the last few years. Conservative schemes, finite volume, 
and artificial density have all contributed to accelerating progress in this 
area. Each has been applied to the analysis of external flow over isolated 
airfoils in two and three dimensions with promising results. 

In the area of internal flows the need to design lighter, more 
efficient turbomachine components has led to consideration of machines with 
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fewer stayes, each with blades capable of higher speeds and higher loading. 
As speeds increaset the numerical problems inherent in the transonic regime 
have to be resolved. These include the calculation of imbedded shock 
discontinuities and the dual nature of the governing equations, which are 
elliptic in the subcritical flow regions but become hyperbolic for 
supersonic zones. The development of supercritical blade design has 
eliminated some of these difficulties by returning to the hoaograph plane 
and explicitly designing blades that will be shock free. However, unlike 
the case of external wing aerodynamics, undesirable off-design conditions 
can develop much faster in turbomachinery flows and are not readily 
controlled or easily corrected. For these reasons the study of off-design 
conditions (high incidence, shocks, etc.) will remain a major requirement 
for Internal flows. In addition, the periodic nature of the geometry for 
turbomachine components adds another degree of complexity which contributes 
to the importance of grid generation as another area of study in its own 
right. 

The present analysis combines some of the most promising transonic 
analysis techniques to calculate tne flowfield surrounding a cascade of 
airfoils. The full potential equation in conservative form has been 
discretized at each point on a body-fitted periodic mesh and a mass balance 
calculated through the finite volume surrounding the point. Each local 
volume is corrected in the third dimension for any change in streamtube 
thickness along the streamtube. Each is also corrected in the tangential 
y(re) direction due to any change in the radial position of the streamtube 
from upstream to downstream. The final result is a tridiagonal matrix 
formulation involving potential corrections at each grid point as the 
unknowns. The right hand side of each equation is equal to the net mass 
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flux through a volume. The resulting system of equations Is solved along 
each grid line by an SLOR procedure. At points where the Mach number 
exceeds unity, the density at the forward (sweeping) edge of the volume Is 
replaced by an artificial density as suggested In reference 1. This 
artificial density Is a function of the real gas censity and Its streamwise 
derivative. 

In the following sections, two body-fitted mesh generators that were 
Incorporated Into tne resulting computer code are discussed. The full 
potential equation In finite difference form Is also reviewed along with 
details on the solution method employed. Subsonic and transonic 
calculations using this method on a NACA0012 profile are compared to other 
computational procedures. Results for several supercritical blade designs 
of Sanz are compared to their hodograpn solutions. Tne quas1-3-U 
capabilities of the method are compared to results from reference 9. 

hesh generation 

The flowfleld Is olscretlzed using a oody-fltteo, periodic 0-type 
mesh. Two techniques were developed to generate computational grios. The 
original grid generation technique, referred to here as the Interpolation 
scheme. Is near-orthogonal and combines the use of analytical functions. 
Interpolations, and simple stretchings. No numerical system must be solved 
for the generation of the grid. The second method, an electrostatic analog 
was developed when the grids generated by the Interpolation scheme were 
found to be unsatisfactory for cascades of high camoer, thickness, and 
solidity. 

Interpolation Scheme 

A summary of this procedure Is shown In figure 1. As suggested In 
reference 2, beginning with the orthogonal, computational plane. 
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-K* < { < +K' 0 < n < K 

the transformation to a physical (x,y) plane as descrioed oy the Jacobian 
elliptic functions sn» cn, and dn(3). Is given as: 


* * 1 ; 
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with K, the complete elliptic Integral of the first kind for m,K(m) ana 
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The specified geometric parameters of chord, stagger angle and blade 
spacing are sufficient to perform the analytical transformations oescrloeo 
by equations (1). After computing K and K' for a blaoe spacing of 
s ■ 2ir and for zero stagger, values of n and t are then selected at 
equal Intervals (or at weighted Intervals If the grid Is to oe dense near 
the body or at the trailing edge). Tables of (x,y) for corresponding (t,n) 
points are generateo using equations (1). In the Cartesian plane this gives 
an orthogonal coordinate system about a cascade of flat plates at zero 
stagger. An example Is shown In figure 2. From this, the conformal 
transformation to a circle Is (4), 

J . X » iy . 

where ( • 6 * u In the circle plane, and the stagger, x. Is zero. 

Equation (2) Is easily Inverted to obtain (a,c) tables If 

X . 0 

s « 

which was the justification for the original calculation of K and K* 
noted above. The parameter R In equation (2) represents the locations of 
the upstream and downstream Infinity points In the conformal plane. It Is a 
function of cascade solidity (c/s) ana stagger angle and can be Interpolateo 
from Information In reference 4 for the specified values of stagger and 
spacing. This same R Is used In calculating (a,c) for the zero stagger 
transform. 
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The resulting circle Is then roteted by applying the exponential terms 
Mith the correct stagger and desired spacing, s, to equation (2). The 
calculated values of 2 *n the Cartesian plane then give an orthogonal grid 
about a flat plate cascade at the desired stagger angle and solidity* 

Figure 3 Is an example of such a mesh. 

The final step In generating this computational mesh Is the adjustment 
of the flat plate grid sc that the Innermost ellipse coincides Mith the 
specified blade shape. This Is accomplished by performing a stralghtforwaro 
Newton Iteration along each line of constant t which Intersects the flat 
plate. The three unknowns to be determlneo at the Intersection of each 
C-llne and the oooy are x, y, and n. The three equations required are 
equations (1) and 

y » F(x) 

where F Is a tabular function giving the specified body shape. 

A successful Iteration yields (x,y)jjQ^^ and npooy «here 


for most blade surfaces, or 


0 < n 


body 


< K 


0 < K < 


’'body 


for the pressure surface of blades whose chordline lies outside the 
blade. The coordinates (x»y)j,pjj^ replace the point on the flat plate 
where the {-line ended. All other grid points on this {-line are shifted by 
means of a simple linear stretching function that reflects the reviseo 
length of the {-line: 




„ ’’body 

K - hi 


^*1 " *pl«te^ * Opiate 


(31 


7 


dnd similarly for each y^. Completed grids are seen In figures 4 ano 5 
for several typical blade shapes. This mesh generation scheme has no 
problem with larger stagger angles (»7b*) nor with high solidities. Cases 
with c/s > 1.45 have oeen attempted with good results. However, If the 
turning angle or camber of the blade Is high, while at the same time, 
stagger and solidity are near the limits mentioned above, the grid lines may 
become crowded on the suction surface. This Is common for turbine blade 
geometries. 

Electrostatic Analog 

An Electrostatic Analog has been oevelopeo to generate body-fitteo 
orthogonal grios tor cascaoes of large turning thickness ano solioity. The 
development of this grid generation scheme begins by considering the 
potential field generated by an infinite spacial array of point charges in 
free space. The density of the charges is assumed to alternate between plus 
and minus one. The mathematical expression for this field is given by 

00 

JT (2 - Zy) ' S l-l)”lh (2 - Zq - ins) (4) 

n»-«e 

where jr(z - 2 q) is the complex potential, z * x ♦ iy, Zq the location 
of the 0th charge, i « ^nd s the distance between charges. 

Equation (4) can be evaluated in closed form by utilizing the theory of 
calculations of residues (ref. 5). Tne result is 

* - (z - z^)» 

^ (-D^lnlz - z^ - ins) « In tanh ^ — (5) 

n»—oo 

(Note that this expression is periodic in y with period 2s.) The 
potential fie d generated by distributing these fundamental solutions over 
the surface of an airfoil in the cascade is given by the integral 
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• f - 2q)02q (6) 

Mhere tl^Q) is the yet unknown source density distribution. On the 
surface of the airfoil (i.e., 2clJ tne real part of m, t* Mill be set equal 
to one. 

Real H(z) « 1 zcL (7) 

At large distances upstream and downstream of the cascade t( 2 ) is assumed 
to approacn zero. This far-field condition imposes the requirement that 


Im y* t(2u)d2j 


(8) 


Equation (6) evaluated on the airfoil surface in conjunction with equations 
(7) and (8) form the integral equation 


1 « Real 


{/ 


y(ZQ)X{Z - 2 Q) 02 rt 7 t 2eL 


*o| • 


0 ■ Im 




IS) 


for the source density function. Once v Is known, the entire potential 
field can in theory be computed by direct integration of equation (b). 

Equation (9) is a singular integral equation, which can be solved by 
paneling procedures similar to those used in solving potential flow 
problems. To employ these procedures one must first factor out of the 
integral equation its singular behavior. For this problem the factor form 
of equation (9) is 
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U - Zq)^ 


♦ # t(2g)1n tanh -~ - g ..- r . - 


( 10 ) 


Next the airfoil surface is divided into a series of lina se9inents. Over 
each of these segments and ln(.2s/(2 - 2o)*3tanhi,(2 - 2 q)«/2sj 

are approximated by their value at the midpoint of the segment. Tne 
singular part of equation (lu) is thus approximated as 
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with the regular part approximated by 
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The auxiliary condition is approximated by 
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(13) I 
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Introducing equations (11), (12), ano (13) Into (10) and restricting the 
value of 2 to the rold-point of each line segment leads to a system of 
linear algebraic equations for the source density olstrloutlon. The 
solution of these equations can be obtained by standard matrix methods, 
with Y known, the complex potential field surrounding the cascade can be 
determined from equation (6). The real part, n and the Imaginary part, t 
of this field form a periodic orthogonal booy-fitted cooroinate system. A 
typical coordinate system Is shown in figure 6. The contours n « constant 
enclose the airfoil, while the curves t« constant project from the 
airfoil to the perioolc oounoary n « U. The line extending to upstream ano 
downstream Infinity are t ■ 1 L. end { « 0, respectively. 

Presently this grid generation scheme utilizes the direct Integration 
procedure outlined above for constructing the ooundlng curvilinear 
cooroinate. for the other grio lines, on the Interior of the mesh, however, 
the direct procedure proved to be Inefficient because of slow convergence. 

To overcome this difficulty an Inverse procedure was developed in which t,n 
become the independent variables and x,y the dependent variables. This 
Inverse procedure is constructed by first noting that the field conversion 
for the electrostatic potential and flux is the Laplace equation 

0, - 0 (14) 



Inverting these equation yields 
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(15) 


as the governing equations for the Inverse mapping (1*e*t x • x(C,n)* 
y • yU»n)). 

The solution of the field equation (IS) suojected to oounoary 
conditions consistent with the coordinates of the oounding curvilinear 
curves yields the desired Interior coordinate curves. A second sample grid 
for a highly comoineo thicx turbine blade generated by the above procedure 
Is shown In figure 7. It Is seen that the grid is hignly regular with 
clustering of grid points In regions of rapio flow variations. The 
periodicity of the grid is quite apparent ano leaos to simple numerical 
approximation of flow variables In the neighborhood of the periodic 
boundary. Flow computations applying this grid generation procedure are 
presented later in the paper. 

solution METhOO 

Governing Equations 

The equation to be solved Is the 2~0 steady, compressible 
full-potential equation in its conservative form 


* (ppy)^ - 0 

Simila<^ to the approach taken by several others (7, 12), this equation Is to 
be written at the finite volume oieroent which surrounds each grid point. 

With superscript n referring to the artificial time step, a typical 


formation at each grid location (1,j) is 
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witn, 

then 

• 0 - U0") ♦ L(««^") . 7(p^”) ♦ -LCi#") (17) 

or 

'•K.j)--»>*‘">t.j--y^ 

where is the volume of the local element, u is the over-relaxation 

parameter, a0!^ j is the unknown change in the local potential field 
af'd j is the residual net mass flow out of the finite volume sur- 
rounoing the grid point. Expanding equation (17) in finite difference form; 



Each volume element is defined by the four surfaces parallel to and locateo 
halfway between neighboring grid lines as shown in figure 8. The tnickness 
of the element (i.e., relative local streamtube height) is obtained from a 
specified height distributiofi or is set to unity for strict two-dimensional 
flow. Variations in the radial (spanwise) position, r, of the streamtube 
are accounted for by scaling all tangential r® coordinates, (y) with a 
ratio of local radius, r, to a given upstreain radius, r . Then, over the 
area, defined by the four surfaces of the volume element 
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Inhere h Is the local stresntube height normalized to unity far upstream. 
The velocities, u, are normal to the element boundaries and located at the 
boundaries. Velocity, II, and density, p, are calculated using the values of 
potential at the previous iteration. 

Using equation (17) and isolating all unknown potential corrections on 
the left-hand side, equation (18) can be expressed in the form 


*^ijn " ^i,j 


i . + A . . , 6^ 1 

i.J i.J-1 i.d-1 




( 20 ) 


where 

and the coefficients 

A « A(p",at,An) (22) 

By writing equation (20) for each point on a grid line of constant t 
(j B 1,J) a tridiagonal matrix system is formed which is then solved by the 
standaro Thomas algorithm (ref. 6) for tne potential corrections 
New values of the field potential are obtained using equation (16; after all 
lines are solved (i 1,1). 

Sweep Strategy 

The grid lines are solved in order beginning with the line which goes 
to downstream infinity and then proceeding clockwise around the blade 
shape. This order results in a solution which sweeps against the flow 
direction on the bottom (pressure) surface of a compressor blade at positive 
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incidence, where the flow is more likely to be subcritical. This minimiies 
the chances of sweeping against a supersonic flowfield. It became evident, 
when symmetric supersonic flows were successfully solved, that direction of 
sweep is not critical. The apparent reason is that additional damping is 
applied by the elliptic-like nature of the artificial density used in 
combination with the sweep strategy. Hore detail on this effect is 
described in the next section. 

Artificial Density 

In equation (18), there is no explicit artificial viscosity term 
necessary for stability in supercritical flow regions and for smearing 
shock jumps. As suggested by Hafez, et al. (ref. 1), an implicit viscosity 
term is incorporated in the density. This modified density is used instead 
of the physical density in equation (18). The form of the modification 
is (1) 

p « p - up^ aS (23) 


where u is the switching function to eliminate the density modification in 
subsonic regions (M < 1) 


V 


max^O, 


1 - 



(24) 


and aS is the total differential component of the density in the flow 
direction, S, defined by the contravariant velocity components a and b. 
Then; 


where. 


iip< 


as 


a 

w- r P an 
a q n 




(26) 


q • ♦ v^ 
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The switching function was split Into two terms so that the effect of each 
component of the density change could be even more efficiently gauged In 
areas of high flow gradient. In addition, an upwind averaged Nach number 
was used: 



*^1+1/2, j±l^ 


I (”in/2,j ^ ”i+l/2±l, j) 


> (2t>) 


Similarly, Upwind differencing Is used for approximating and 
0 ^. The physical densities used to calculate and p^ are updated 
as each line Is reacheo. Implicit In this arrangement Is the creation of a 
timelike cross derivative of density, p^^, when sweeping against the flow 
direction. This Is because tne density at the current line has been 
calculated using the values of potential and velocity at the latest 
Iteration (1,e., artificial time step), while the values at the next 
upstream line, where the upwind differences are formed, remain at the 
previous time step until U In turn Is solved. However, when sweeping with 
the flow, the upstream line from which the density gradients are formed has 
already been solved and updated. Thus this additional damping due to the 
time lag of densities Is not present when solving with the flow and the 
artificial compressibility reduces to the original space-lag only. 

Boundary Conditions 

The boundary conditions applleo to the solution system Include a 
specified upstream axial mass flow in the form of an oncoming Mach number, 
plus Incidence and stagger angles. Using global continuity. 


( 2 /) 
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where A Is the flow erea nonna) to the streamtube. 

The specified flow angle at downstream Infinity together with a simple 
Iteration procedure on the Isentropic density ratio gives the downstream 
Mach number boundary condition. For a stagger angle, x, and upstream 
Incidence (chordwise) of equation (27) can be rearranged to yield 



where 


a • o ♦ X 

Using a straightforwaro Newton iteration on equation (2b) yields M^. From 
this tne relative soeed ratio is available as: 


q M a 

90 




(29) 


Periodicity is maintained in the tangential (y) direction by enforcing the 
relation 

d(x,y ♦ s) « d(x,y) + C (30) 

where s is the blade spacing and 

y ^(x,y+s) 

) V . dl . V s (31) 

(x,y) 

On the grid line defining the body, a standard zero normal flow 
condition is imposed 


V# • n - 0 


(32) 
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RESULTS 

Some results obtained from these methods are shown In figures 9 
through 21. Verification of tne code at low speed Is provided In figure 9 
by the comparison to the methods of Caughey (ref. 7) for the NACA0012 blaoe 
shape. The transonic nonllfting case Is shown In figure 10 for this same 
blade. The shock strength and location agree well with similar results 
using the methods of reference 7. Tne solution Is surprisingly symmetric, 
particularly in light of tne fact that tne solution on the lower surface was 
obtained by marching against the supersonic flow and the upper surface was 
obtained by marching with the flow. This relative shift In marching 
direction occurs because the solver was written such that it makes one 
continuous clockwise sweep of all grid lines around the airfoil. Thus there 
Is a primary difference between this method and other transonic tecnniques. 
It is apparently insensitive to the relative direction of the supersonic 
flow. For the cases tested thus far, this insensitivity to flow olrection 
does not seem to oe affected oy mesh size (down to aX » 2 percent chord). 

The symmetry of the solution can be seen very effectively in the isomacn map 
of figure 11. 

The sensitivity of the system to a near-shock condition car. be seen In 
the results for the supercritical blade shown In figure 12. This Is a 
compressor stator tip section designed for NASA-Lewis by Sanz using a method 
based on Bauer, et a1. (ref. 8). The trailing edge ends In a cusp. The 
fast recompression is captured well without any over-reaction or steepening 
Into a shock. 

The capability to approximate moderate 3-0 streamtube effects Is 
demonstrated In figure 13 wnere both the strict 2-0 flow and quas1-3-0 flow, 
as corrected for a given streamtube convergence profile, are compared for 
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the double circular arc stator of figure 5. Agreement is excellent with 
results using the method of reference 9. For this case there was a 
7-percent gradual contraction in the streamtube from far upstream to far 
downstream. In addition, the streamtube varied its position along the span 
(radial shift) by 0.7 percent over the same distance, a negligible effect. 

The resulting isomach map is shown in figure 14. 

A comparison of results for a typical turbine blade section is shown in 
figure 15. The blade used is designated as the AACE(II) airfoil by 
R. A. Delaney, to whom we are grateful for supplying the cooroinates and 
experimental data (ref. 10). Agreement is good on tne pressure side, but 
the calculated envelope is slightly low on the suction surface. This is 
consistent with the metal blade coordinates used in this inviscid 
calculation. If a suitable boundary layer analysis were incluaeo and 
displacement thickness added to the profile, the suction surface calculation 
would more accurately matcn the experimental data. The location of each 
peak in speed is predicted accurately, even near the trailing edge, despite 
the inherent difficulty for potential flow methods to go around sharp bends 
ct high velocity without some artificial limitation on calculated speeds. 

Such limits were not necessary for the coarse mesh used for figure 15 (26 
points around the blade). However, finer meshes were unsuccessful without 
some artificial barrier imposed on Mach number. 

Figure 15 also gives an interesting comparison between the two mesh 
generators described earlier. The results for both grids are virtually 
identical in the subsonic region as would be expected. In the high-speed 
region, the exact-orthogonal mesh (electrostatic analog) appears to be 
slightly more conservative overall than the near-orthogonal grid. 

Subsequent experience with ether denser meshes showed that tne 

■I 

. j 
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exact-orthogonal mesh provided more stability In high gradient areas, such 
as the trailing edge. Figures 16 and 17 give views of the two grids used 
for the calculations of figure 15. An Isomach map for the AACE(II) turbine, 
using the near-orthogonal grid. Is shown In figure 18. 

As a final example of tne Importance of streamtube effects on 
supercritical blade design, figure 19 shows a series of calculations 
performed on a thick compressor stator hub section developed for NASA-Lewls 
by Sanz (via ref. 8). The present method was first used to calculate the 
strict 2-0 flow over the blade. This is seen to agree very well with the 
Sanz hodograph design envelope except at the trailing edge, where the 
Idealized blade of Infinite length (hodograph solution) was replaced b> a 
cap of constant radius to connect pressure and suction surfaces for the 
present method. The remaining two curves are quasi-3-D cases which 
hightlight tne effect of changes in spanwise (radial) position of the 
streamtube. Both cases involve a streamtube thickness cnange of about 
-14 percent (AVDR » 1.15). The axial distribution of this thicxness change 
(h/h_^ as well as tne distribution of spanwise streamtube position 
(r/r_j as shown in figure 19 were obtaineo for this Sanz hub section by 
using the J(at>anis meridional plane computer^rogram (ref. 11). The results 
are typical for a section near the hub wall with thickness decreasing, and 
spanwise position Increasing about 5 percent. One quas1-3-0 case In 
figure 19 has only the thickness change applied, with radius held constant. 
The other case nas both thickness and radius cnange. Note that the 
Increasing radial position of the streamtube has a strong decelerating 
effect on the flow since the blade-to-blade passage now diverges (re larger 
at trailing edge than at leading edge). A oesign using this result could 
tolerate higher blade flow rates than one which ignored these radial effects. 
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Figure 20 shows the Isomach map for the case of thickness variation 
only. The shock appears as expected near 40 percent chord. Figure 21 Is 
for the case of both thickness and radius variation. Note that the shock 
has washed out and flow In the passage has decelerated. 

CONCLUSIONS 

A reliable method has been presented for calculating the flowfleld 
about a cascade of arbitrary 2-0 airfoils which can be made to approximate 
the three-dimensional flow in a turbomachinery blade row oy correcting for 
streamtube convergence and radius shift. The method Is a fully conservative 
solution of the full potential flow equation using a body-fitted periodic 
mesh, with an artificial oensiiy imposed in the transonic region to insure 
the capture of shock waves. A computer code has been developed using the 
methods described aoove. It has been tested using a variety of 
turbomachinery and general blade profiles. The calculation of a shocked 
symmetric NACA blade has verified that this method is insensitive to the 
relative Directions of supersonic flow and SLOR sweep. Comparison of 
results for several well-known supercritical blades has shown satisfactory 
agreement with their hodograph solutions. The approximation of moderate 3-0 
streamtube effects was demonstrated by a good comparison between this method 
and results using the Katsanis code on a OCA compressor stator. 2-0 and 
quasi-3-0 calculations for a supercritical compressor stator Indicated that 
the moderate 3-0 effects included in this method are important enough to 
radically affect design considerations in high-speed regions of the blaoe. 
These results all Indicate that the method described is a useful one for 
analyzing off-design performance of turbomachinery blading. 
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APPENDIX A 
LIST OF SYMBOLS 

(eqs. (^0)(22)} vdriable coefficient of change In potential for 
field equation 

(eqs. (27) (28)) flow area normal to the streamchannel direction 
(eq. (19)) total surface area of finite volume element 
contravariant velocity component (n direction) in computational 
space 

contravariant velocity component (t direction) in computational space 
circulation about tne airfoil as defined by the given values of 
flow angle 
blade chord 

length of passage along x-axis, from leaoing edge to trailing edge 
when blade is at setting angle 
tabular function of body shape (x,y) pairs 
transposed terms in field equation involving known (previously 
determined) potential corrections 
local streamtube height normalized to unity far upstream 
complete elliptic integral of the first kind for (k^) 
dt/Vl - k^ sin^t 

complete elliptic integral of the first kind for (1 - k^) 
complex electro-static potential field generated by infinite array 
of point Charges 

g-(»rc)/(s cos x) parameter, defined by cascaoe geometry, for 
obtaining Jacobian elliptic functions 
full potential operator (eq. (17)) 

length of bounding surfaces for finite volume (see fig. 8) 
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’ M Mach number 

m 

n unit vector, normal to airfoil surface 

q magnitude of flow speed, normalized to unity at far upstream 

* 

conditions 

‘ R residual net mass flow out of a finite volume using potential 

values from previous iteration (eq. (ly)) 

' St R G, total residual right hand side for field equation at each 

volume 

r local spanwise radial position of streamchannel normalized to unity 

far upstream 

S flow oirection as defined by contravariant velocities a and b in 

I equation {2S) [] 

s blade spacings 

T a combination of Jacooian elliptic functions, defined in equa- 

tion (1) 

u x-componer.t of velc:ity in the physical plane 

A A 

V total physical velocity vector, ui + vj (eq. (31)) 

V y^component of velocity in the physical plane 

X streamwise position coordinate, i.e., in figure 4 x-axis Is 

hori zonal 

y tangential position coordinate, i.e., in figure 4 y-axis is 

vertical 

z coordinate to transform a physical point to complex plane, 

z « X iy 

2^ cooroinate of beginning of one line segment. Blade is divioed into 

segments in solving for electrostatic analog grid 
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2|j coordinate at the end of a segment of the blade surface 

Zq coordinate of zeroth charge for electrostatic analog 

o Incidence angle of flow as measured from chordllne of blade 

0 angle of flow as measured from x-axIs 

t ratio of specific heats 

y(Zq) (electrostatic analog grid procedure). Source density distribution 
r defined in equation (1) 

6 abscissa coordinate In the plane of tne conformal circle transfor- 

mation (Interpolation scheme, eq. (2/) 
aP change In the local potential field, the unknowns In equation (18) 

c ordinate In conformal transformation plane (eq. (2)) 

t coordinate in computational space whose contours In physical space 

radiate from the blade and terminate on the periodic boundary 
n coordinate in computational space whose physical contours surround 

the airfoil 

« tangential coordinate in radians, y * re 

X stagger or setting angle of the blade row, as measured from x-axis 

to chord line 

II switch function to control the artificial density modification 

(eq. (24)) 

« 3.14159265 

p physical density 

artificial density as described in equation (23) 

P flow velocity potential function 

H electrostatic potential function due to cnarge distribution on 

blaoe surface (eq. (6)) 

w over-relaxation parameter, 1 < w < 2 

Q defined In equation (1) 
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Subscripts: 

1 position on computational plane In {-direction 

j position on computational plane In ^-direction 

X differentiation with respect to streamwlse position coordinate, x 

> olfferentlatlon with respect to tangential coordinate, > 

{ differentiation with respect to { 

1 ) differentiation with respect to n 

-00 conditions at far upstream station 

♦oo fjr downstream conditions 

Superscripts: 

n artificial tim- step (iteration level) 

-(bar) (eq. (19) and fig. 8) average over the length, t, of the velocity 
component normal to the element surface along t 
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Figure 1(a). - Analytical transformation steps for near-orthogonal grid development 



body shape 

Figure Kb). - Interpolation and shear transformation steps for near-orthogonal grid development 
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Figure 4. - Near-orthogonal mesh for a symmetric blade. 



Figure 6. - Typical mesh generated by electro-static analog procedure 
fora turbine blade with trailing extension. 



Figure 7. - Computational mesh developed by electrostatic analog 
procedure for another thick turbine blade. 








supercritical stator tip section 









Figure 15. - Comparison of experimental and calculated 
results for turbine cascade, AACE (III. 
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Sanz. hodograph design curve 
Present method, strict 2D 
Present method, ho^h^- ,849, r^r. 
Present method, hoo/h^* .849, roo/r 
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jure 19. * Effect of streamtube’s radial position and thickness on flow about a 
supercritical compressor stator hub section. 








